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Abstract 

Spline Galerkin approximation methods for the Sherman-Lauricella inte¬ 
gral equation on simple closed piecewise smooth contours are studied, and 
necessary and sufficient conditions for their stability are obtained. It is shown 
that the method under consideration is stable if and only if certain operators 
associated with the corner points of the contour are invertible. Numerical 
experiments demonstrate a good convergence of the spline Galerkin methods 
and validate theoretical results. Moreover, it is shown that if all corners of 
the contour have opening angles located in interval (O.lvr, l.Ovr), then the cor¬ 
responding Galerkin method based on splines of order 0, 1 and 2 is always 
stable. These results are in strong contrast with the behaviour of the Nystrom 
method, which has a number of instability angles in the interval mentioned. 
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1 Introduction 


Let D be a simply connected planar domain bounded by a piecewise smooth curve 
r. It is well known that the solution of various boundary value problems for the 
biharmonic equation 


A'^u{x,y) 


d^u d^u d^u 

dx^ dx'^dy'^ dy^ 


0, {x,y)eD, 


where A is the Laplace operator, can be constructed via solutions of boundary 
integral equations. Consider the biharmonic Dirichlet problem 

A^|d = 0, 


w|r — /i, 


du 

dn 


/2, 


( 1 . 1 ) 


where du/dn denotes the normal derivatives and /i, /2 are sufficiently smooth func¬ 
tions dehned on the boundary L. Setting z = x + iyO"^ = —1, one can identify 
D with a domain in the complex plane C. This problem arises in various appli¬ 
cations, in particular while considering the behaviour of viscous flows with small 
Reynolds numbers, bacteria movement, deflection of plates, elastic equilibrium of 
solids, sintering P [TTl [T71 [20l |2^ [23l l2l] . 

Let us equip the curve L with the counterclockwise orientation and consider the 
Sherman-Lauricella equation 


‘"W + i //K) (H) ■ ^ (H) ^ ^ 

( 1 . 2 ) 

where the bar denotes the complex conjugation and oo is an unknown function. 
Equation fll.2p originated in works of G. Lauricella (see [H]). He was the hrst who 
used the method of integral equations in elasticity. Later D.I. Sherman rewrites 
Lauricella equation in a complex form and proposes a new simple way to derive 
it [27]. The equation fll.2p is uniquely solvable in appropriate functional spaces, 
provided / satishes certain smoothness conditions and 


Re J f{t)dt = 0, (1.3) 

[131 [20l [21]. Moreover, let a = a{x,y), {x,y) G T denote the angle between the 
real axis M and the outward normal n to T at the point {x, y) and let 1 be the unit 
vector such that the angle between 1 and the real axis is a — 7r/2. If one dehnes the 
function / = f{t) = f{x, y), t = x + iy by 

(/ 2 (f) + ^ ^ r. 


fit) ■= e ' 


(1.4) 
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then the solution of the Sherman-Lauricella equation fll.2p with such right-hand 
side / can be used to determine a solution of the boundary value problem fll.ll) . 
More precisely, if w is a solution of the equation fll.2p with the right-hand side fll.dp . 
consider two holomorphic functions (p = (p(z) and ^|J = z E D dehned by 


ip{z) = 

ip{z) = 

According 

problem 

Therefore, 






f (‘^(0 


27ii Jr C ~ z ^ 27ri Jr C ~ z 27ii Jr (C — z)"^ 


(1.5) 


d(, z E D. (1.6) 


to [20], the pair {(p{z),'il’{z)} represents a solution of the boundary value 

(p(t) + tpi'it) + '0(t) = e“*“ ^ e r. 

by [m Lemma 5.1.4] the function 

u{x,y) := Re {z(p{z) +'ip{z)), z = x + iy E D (1.7) 


is the solution of the boundary value problem fll.ip . 

Thus if an exact or an approximate solution of the integral equation fll.2l) is 
known, a solution of the biharmonic problem fll.ip can be obtained by using for¬ 
mulas fll.5p . fll.Op and fll.7p . Therefore, the main effort should be directed to the 
determination of solutions of the Sherman-Lauricella equation fll.2p . Note that the 
Nystrom method for the Sherman-Lauricella equation on smooth contours has been 
used in [HJ [T8| to hnd approximate solution of biharmonic problems arising in fluid 
dynamics. However, the authors of those works have not presented any stability 
conditions for the method considered. If T has corner points, the stability study 
becomes more involved since the integral operators in fll.2p are not compact. For 
piecewise smooth contours, conditions of the stability of the Nystrom method are 
established in PIZ]. These results have been used in [8] in order to construct a very 
accurate numerical method to hnd solutions of the biharmonic problem fll.ll) in 
piecewise smooth domains in the case of piecewise continuous boundary conditions. 

In the present paper, we consider spline based Galerkin methods for the equation 
fll.2p and study their stability. It is shown that the corresponding method is stable if 
and only if certain operators R'^ from an algebra of Toeplitz operators are invertible. 
These operators depend on the spline space used and on the opening angles of the 
corner points r G T. Unfortunately, nowadays there is no analytic tool to verify 
whether the operators in question are invertible or not. Nevertheless, we propose a 
numerical approach which can handle this problem. Thus spline Galerkin methods 
are applied to the Sherman-Lauricella equation on simple model curves and the 
behaviour of the corresponding approximation operators provide an information 
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about the invertibility of the operators R^, r G F. Note that in comparison to 
the Nystrom method, the implementation of spline Galerkin methods to solve the 
Sherman-Lauricella equation, requires more preparatory work. On the other hand, 
numerical experiments suggest that these methods have no ’’critical” angles located 
in the interval [O.lvr, l.Qvr], i.e. if the boundary F does not possess corners with 
opening angles from the interval mentioned, then these methods are stable. In a 
sense, this is similar to the behaviour of the corresponding approximation methods 
for Sherman-Lauricella and Muskhelishvili equations in the case of smooth curves 
which always converge laiiiCQ]. Of course, one also has to study the opening angles 
in the intervals (0,0 .Itt) and to (l.Ovr, 27r) but this is a time consuming operation 
and will be considered elsewhere. 

2 Splines and Galerkin method 

We start this section with the construction of spline spaces on the contour F. Let 
7 = 7 ( 5 ), s G M be a 1-periodic parametrization of F, and let A4r denote the set of 
all corner points tq, ti, ..., Tq-i of F. Without loss of generality we can assume that 
Tj = ■yQ/q) for all j = 0,1,..., g — 1. In addition, we also suppose that the function 
7 is two times continuously differentiable on each interval {j/q, {j + l)/q) and 



Note that the last condition is not very restrictive and can always be satisfied by 
changing the parametrization of F in an appropriate way. 

Let / and g be functions defined on the real line M, and let f * g denote the 
convolution 



of / and g. If y is the characteristic function of the interval [0,1), 



1 if s G [0,1), 
0 otherwise. 


then (j) = refers to the function defined by 



Recall that for any given non-negative integer d, the function 0 generates spline 
spaces on M. Thus if an n G N is fixed, then closure in the L^-norm of the set of all 
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finite linear combinations of the fnnctions 4>nj{s) '■= 4>{ns — j), j E constitntes a 
spline space on M. 

Using the above dehned spline fnnctions, one can introduce spline spaces on the 
contour F. More precisely, for a hxed non-negative integer d and an n G M, n > d+1, 
we denote by = S^(T) the set of all linear combinations of the functions 

^njit) :=^{ns - j), f = 7 (s) e r, j = 0,1, ■ ■ ■ ,n - {d + 1), s G M, 


the support of which belongs entirely to one of the arcs [r^, r^+i), k = 0,... ,q and 
Tq+i := tq. This dehnition is correct since the support supp0 of the function 0 is 
contained in the interval [ 0 , d -|- 1 ] [26] and 7 is a 1 -periodic function. 

In what follows, we also consider operators acting on various subspaces of the 
Hilbert space U = of all sequences {^k) of complex numbers ^k,k eZ satisfying 
the condition 


llfe)l|: = 



< 00. 


The space U is closely connected to spline spaces on the real line M. Thus the 
following result is true. 


Lemma 2.1 (0]) Let n G N. Then there are constants Ci and C 2 such that for any 
sequence {fQ E the relations 


||(^fc)|| < ci^/n\ T 




L2 (k) 


< 


C2 


n 


11 ( 6)11 


hold. 


Further, let L^{T) denote set of all Lebesgue measurable functions / such that 

II/IIl 2 := (^ |/(t)|^ds) <CX), 

and let Hr : T^(r) -E T^(r) be the operator corresponding to the Sherman-Lauricella 
equation fll.2p . It is well known that the operator A-p is not invertible on the space 
L^(r) [20]. On the other hand, the invertibility of the corresponding operator is a 
necessary condition for the applicability of any Galerkin method to any operator 
equation. Therefore, for approximate solution of the equation 01 . 21 ) we use the 
equation with the operator Bp instead of Ap and chose the right-hand sides / of the 
initial equation fll.2p from a suitable subspace of L‘^{T). More precisely, let lU^’^(r) 
denote the closure of the set of all functions / with bounded derivatives in the norm 


lyi 


\f{t)\^ds+ / \f{t)\^ds 


1/2 
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and let Tr : refer to the operator defined by 


Truj{t) : = 


{t — a) 2m Jr \{C ~ «) 




(C - ay 


(2.1) 


where a is a point in D. 

Theorem 2.1 (|I7]) IfT is a simple closed piecewise smooth contour, then operator 

Br ■= Ar + Tp 

is invertible on the space L 2 (T). Moreover, if function f G Wf{T) satisfies the 
condition ca, then the solution of the equation 


BrUJ = f 


( 2 . 2 ) 


belongs to the space Wf{V) and is a solution of the original Sherman-Lauricella 
equation (O. 

Thus if the right hand sides / G Wf{T), the corrected Sherman-Lauricella equa¬ 
tion can be used in order to hud an exact or an approximate solution of the equation 
fll.2p . In the present paper, we employ spline based Galerkin methods to the equa¬ 
tion fl2.2l) and study their stability and convergence. Let us describe these methods 
in more detail. First of all, we normalize all the basis spline functions used. If n is 
hxed, then for any j G Z the norm ||(/>nj|| of any basis element (fnj is 


110: 


nj I 


= — [ ^(s) ds. 

n In 


Therefore, if refers to the number 


■= 


^(s) ds. 


then 


(2.3) 

(2.4) 


0nj • \/^ 0nj) J ^ ^ 

are unit norm vectors. An approximate solution of the equation fl2.2p is sought in 
the form 

^nijJ ^ ^ ®fc0nfc(0) (^•^) 

Kk&si{v) 

the coefficients of which are obtained from the following system of algebraic 
equations 

(5rWn,0ni) = (/,0nj), 0nje^^(r). (2.6) 
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An important problem now is to stndy the solvability of the equations fl2.6l) and con¬ 
vergence of the approximate solutions to an exact solution of the original Sherman- 
Lauricella equation fll.2p . In Section [3l this problem is discussed in a more detail 
but, at the moment, we would like to illustrate the method under consideration by 
a few numerical examples. Thus we present Galerkin solutions of the equation fll.2p 
with the right-hand side / = /i, 

fi{z) = f{x,y) = - 12xy‘^+ i{4y^ - Ux'^y); z = x + iyeT, (2.7) 

on the unite square and rhombuses, and trace the evolution of the solution when 
the initial contour is transformed from the unit square into rhombuses with various 
opening angle a. Some of these contours have been used in in order to illustrate 
the behaviour of the Nystrom method. Note that in the corresponding examples 
from [7j, approximate solutions of the equation fll.2p with the right-hand side 

f2{z) = 1^1 

have been determined. We apply the spline Galerkin method to the equations with 
such right-hand side, too. The results obtained have a very good correlation with [7] 
and the error evaluation for both cases are reported in Table [H where El\ denotes 
the relative error ||a; 2 n—i^n|| 2 /||i^ 2 n ||2 computed for the righthand side /* and equation 
fll.2p is considered on the rhombus with the opening angle a. In addition. Figures 


Table 1: Relative error of the spline Galerkin methods 


n 

/9 

71,71/2 

71,77/6 

E^^ u 

71,77/4: 

71,77/0 

E^" /o 

71,77/2 

71,77/6 

71,77/0 

128 

0.0373 

0.6194 

1.3577 

2.1716 

0.0121 

0.0217 

0.0205 

256 

0.0198 

0.0268 

0.2046 

0.6169 

0.0067 

0.0112 

0.0245 

512 

0.0096 

0.0059 

0.0616 

0.1888 

0.0045 

0.0102 

0.0193 


[IHllshow the convergence of the approximate solutions of the equation fll.2p with the 
right-hand side fl2.7p obtained by the Galerkin method based on the splines of degree 
d = 0 and the transformation of these approximate solutions when n increases. 

Let us mention a few technical details related to the examples below. Thus the 
rhombus with an opening angle a is parameterized as follows. 


7(s) = < 


V 


4s — cos 

(4s - l)e*“ -isin 
-(4s-2) + cos(^|) 
— (4s — 3)e*“ -f isin ^ 


gia/2 


gia/2 


^ia/2 


if 0 < s < 1/4, 
if 1/4 < s < 1/2, 
if 1/2 < s < 3/4, 
if 3/4 < s < 1. 


a 

2 


( 2 . 8 ) 

























GALERKIN METHODS FOR SHERMAN-LAURICELLA EQUATION 


Moreover, we have to compute the scalar products {BrOJn, 4>nj)- Recall that supp (pnj c| 
{j + d + l)/n] and use the Gauss-Legendre quadrature rule with quadrature 
points which coincide with the zeros of the Legendre polynomial P 2 a{x) on the canon¬ 
ical interval [—1,1], scaled and shifted to the interval [j/n, {j + d + l)/n]. More 
specihcally, the corresponding formula is 


rij+d+l)/n _ 24 

{BrUJn,4>nj) = / BrUJn{(l>{s))4>nj{Hs))ds ^ y^WfcRr^n(0(gfc))0nj(0(gfc)), 

k=l 

(2.9) 

where Wk, Sk are the Gauss-Legendre weights and the Gauss-Legendre points on the 
interval [j/n, {j + d + l)/n]. In order to hnd the values of the corresponding line 
integrals at the Gauss-Legendre points, the composite Gauss-Legendre quadrature 
is used O Section 3], namely. 


k{t,T)x{T)dT = / k{'y{a),'y{s))x{'y{s))'y'{s)ds 


m—1 r—1 

Wpk{-f{a),-f{sip))x{-f{sip))Tlp/m, 

1=0 p =0 


( 2 . 10 ) 






Figure 1: Approximate solution uJn(t) of the Sherman-Lauricella equation fll.2p on the 
unit square F with / := fi defined by 02.71) and <7 = 0. From the left to the right: 

n = 128,256,512,1024 
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where = 'y'{sip) with m = 40 and r = 24. 

Table [Hand Figures HUH show a good convergence of approximate solutions if the 
corner point of the contour has an opening angle close or equal to 7r/2. On the other 
hand, the presence of opening angles of a small magnitude can cause problems and 
lead to a convergence slowdown (see Figures [31ll]). Note that although the focus of 
this work is on the stability, the error estimates presented in Tabled] are comparable 
with estimates of recent work [16] for fast Fourier-Galerkin method for an integral 
equation used to solve boundary value problem fll.ip in smooth domains. Moreover, 
further improvement of the convergence rate is possible if for the approximations of 
singular integrals and inner products arising in the Galerkin method one employs 
graded meshes of various kind urn- 


3 Galerkin method. Local operators and stability 

Our next task is to hnd conditions of applicability of the spline Galerkin methods to 
the equation fl2.2l) . It is worth mentioning that for smooth contours F, the methods 
considered here are always applicable and provide satisfactory results. For details 
the reader can consult [5] , where similar methods for the Muskhelishvili equation on 
smooth contours are considered. On the other hand, the presence of corners changes 
the situation drastically, and the applicability of the approximation method is not 
always guaranteed. 

Let Pn be the orthogonal projection from T^(F) on the subspace 5'j((F). Then 
the systems fl2.6p is equivalent to the following operator equations 

PnBrPnl^n = Pnf, n EN. (3.1) 

Definition 3.1 We say that the sequence (PnE^Pn) is stable if there is an m eM. 
and an uq E N such that for all n > uq the operators PnE^Pn : S^{T) Sf^{T) are 
invertible and 

\\{PnBrPn)~^Pn\ \ < m 

for all n > uq. 

Recall that if the stability of the corresponding sequence {PnB-pPn) is established, 
then the convergence of the Galerkin method and error estimates can be obtained 
from well known results, cf. m Section 1.6, inequality (1.30)]. Therefore, in this 
work we mainly deal with the stability and our approach is based on G*-algebra 
methods often used in operator theory. Let CaddiL"^(X)) refer to the real G*-algebra 
of all additive continuous operators on the space L^(F). One can show [H] that any 
operator A E Cadd{L'^(J')) admits the unique representation A = Ai + A 2 M where 
AuA2 are linear operators and M is the operator of complex conjugation. This 
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Figure 2: Approximate solution of the Sherman-Lauricella equation fll.2p on the 

rhombus F, a = vr/S with / := /i defined by fl2.7l) and d = 0. From the left to the 
right: n = 128,256,512,1024 






Figure 3: Approximate solution Unit) of the Sherman-Lauricella equation fll.2p on the 
rhombus F, « = 7r/4 with / := /i defined by (12.71) and d = 0. From the left to the 
right: n = 128,256,512,1024 
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Figure 4: Approximate solution of the Sherman-Lauricella equation fll.2p on the 
rhombus F, a = vr/S with / := /i defined by fl2.7l) and d = 0. From the left to the 
right: n = 128,256,512,1024 


representation allows one to introduce the operation of involution on CaddiL"^(X)) as 
follows 

A*:=Al + MAX (3.2) 

with A\,A2 being usual adjoint operators to the linear operators 741,^42, cf. m 
Theorem 1.3.8 and Example 1.3.9]. By we denote the set of all bounded sequences 
{An) of bounded additive operators A^ : imP„ —> imP„ such that there is an 
operator A G Cadd{L‘^{^)) with the property 

S - limAnPn = A, S - lim {AnPn)*Pn = A*, 

where s — Van. An denotes the strong limit of the operator sequence {An). 

Provided with natural operations of addition, multiplication, multiplication by 
scalars A G C, with an involution introduced according to fl3.2p . and with the norm 

||(Al„)|| := sup llAlnll, 

riGN 

the set becomes a real C'*-algebra. Consider also the subset C A^ consisting 
of all sequences (J„) of operators Jn : imP„ —)■ imPji which can be represented in 
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the form 

Jn = PnTPn + Ti G M, 

where the operator T belongs to the ideal ]Cadd{L^{^)) C Cadd{L'^(J')) of all compact 
operators and the sequence (Cn) tends to zero uniformly, i.e. 

lim ||C„|| = 0. 

n—)-oo 

The stability of sequences from the algebra can be characterized as follows. 

Theorem 3.1 (cf. |11|, Proposition 1.6.3]) A sequence {AA) G .4,^ such that 
A ■= s — hm4„P„ is stable if and only if the operator A is invertible in Cadd{.L‘^{^)) 
and the coset {An) + is invertible in the quotient algebra AF/ . 

Consider now the sequence {PnBvPn) of the Galerkin operators dehned by the 
projection operators Pn- Recall that on the space T^(r) the sequence of the orthog¬ 
onal projections (P„) strongly converges to the identity operator I and P* = Pn, n G 
N. It implies that for any operator A G £add(T^(r)) the following relations 

s - lim PnAPn = 4, s - lim(P„4P„)*P„ = 4* 

hold [25]. Therefore, combining Theorem 12.11 and Theorem 13.11 one obtains the 
following result. 

Corollary 3.1 Let T be a simple closed piecewise smooth curve. The spline Galer¬ 
kin method fl3.ip is stable if and only if the coset {PnBvPn) + is invertible in the 
quotient algebra Af^ / . 

Thus in order to establish the stability of the Galerkin method, one has to study 
the invertibility of the coset {PnBvPn) + in the algebra AF j . This problem 
can be tackled more efficiently, if we restrict ourselves to a smaller algebra containing 
the coset {PnBvPn) + J■ More precisely, let M refer to the operator of the complex 
conjugation, 

M(f){t) := 0(t), t G T, 
and let Sv be the Cauchy singular integral operator, 

Srm ■=- fpF 

Til JvC-t 

Consider the smallest closed real C'*-subalgebra of the algebra A^ which con¬ 
tains all operator sequences of the form {PnMPn), {PnSvPn) and also the sequences 
{PnfPn),f G CR(r) and (G„), where hm„_^oo ||C„|| = 0 and CK(r) is the set of all 
continuous real-valued functions on the contour T. 
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Remark 3.1 It follows from flU ITE. \2^ that C and that the sequence 
(PnBrPn) belongs to BA. Therefore, BA/ is a real C*-subalgebra of JA j , and 
by m Corollary l.f-lO] the coset (PnBrPn) + is invertible in JA/if and 
only if it is invertible in BA/ . Therefore, one can now study the invertibility of 
the coset {PnBYPn) + in the smaller algebra B^/. To this end we will employ 
a localizing principle. 

Thus with each point r G T we associate a model contour T,- as follows. Let Or 
be the angle between the right and the left semi-tangents to T at the point r, and 
let At refer to the angle between the right semi-tangent to T and the real line M. 
Consider now the curve 

:= eAhT+oA]^- u 

where and denote the positive semi-axis M’*' correspondingly directed to 
and out of the origin. Further, on each such contour F,-, r G F we consider the 
corresponding Sherman-Lauricella operator 

Ar = I + Lr- KrM, (3.3) 


where 


in ^ (^) ■ 

Analogously to the algebra B^ and to the ideal one can introduce algebras B^'^ 
and ideals C B^'", r G F, which allow to establish conditions of the applicability 
of the corresponding Galerkin method for the operators (13.3p . For this we also 
need appropriate spline spaces on both the contour F,- and the positive semi-axis 
:= These spline spaces can be constructed by using the functions fl2.4p 
again. More precisely, consider the functions 



Anjis) if t = e'-^^s . ^ 

0 otherwise ^ ~ 

An,j-d{s) if t = eAh^+^As 

0 otherwise 


j < 0 


(3.4) 


Let 5'(((F.r) and S'(((M’'') be, respectively, the smallest closed subspaces of L 2 {Tr) 
and L^(M+) which contains all functions fl3.4p and all functions Anj, J > 0 of fl3.4p 
for At = 0. Moreover, let Pf, n G N and denote the orthogonal projection 
onto subspaces Sf^iTr) and S')((M+), respectively. In order to study the stability of 
the sequence {PfArPf), one can apply Theorem 13.11 and Remark [3.II to obtain the 
following result. 
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Corollary 3.2 The sequence {P^At-P^) G is stable if and only if the operator 
Ar is invertible in and the coset (PfArPf) + is invertible in the quotient 
algebra 

Further, let L|(M+) be the space of all pairs {(pi, <^ 2 )'^, V^i, T 2 ^ provided 

with the norm 

\\iTi,T 2 f\\ ■= (I i</^iir+ 

and let rj : —)■ LKM^) be the mapping dehned by 

r]{ip) = s G 

where denotes the transposition of the vector a. It is clear that r; is a linear 
isometry from onto LKM"*"). Moreover, the mapping 4/ ; Cadd{L‘^(J'r)) 

CaddiLlO^'^)) dehned by 

^{A) = riAri~^, (3.5) 

is an isometric algebra isomorphism. In particular, straightforward calculations show 
that 


vb(F;) = diag(P+,P+), 


where 


d'(M) = diag(M,M), 


= 

= 


0 K. 

Afe. 0 

0 

_pi2{iST+eT)j^ 


0 


= oTT- 


1 1 /•°° 

22m Jq — ((j/s)e*^^ 1 — [a/s)ed‘^'^~^A 


ds 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


Mg^^p{a) 



sin Or 

(1 — {a / s)e^^'^Y 



and the symbol M in the right-hand side of (13.7p refers to the operator of the complex 
conjugation on the space L^(M’''). Moreover, one can observe that the operators Afg^ 
and AA.e^ have a special form - viz. 



K^p{a) : 


s 


(3.10) 
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■■= 'nefiu) 

1 u sin Oj. 

27r 1 — P ’ 

if 

K 

— Afe.,, 

(3.11) 


■= mefiu) 

1 u sin 9r 
vr (1 — ’ 

if 

K 

= Me^. 

(3.12) 


On the space R of the sequences of complex numbers ^k, k = 0,1, ■ ■ ■, 

OO 

P :={(&)£» 

k=0 


the function dehnes a bounded linear operator A{kQQ) with the matrix represen¬ 
tation 






0(m) 


du 

u + q 



1 


where Vd is the constant fl2.3p . 


Theorem 3.2 Let and be the functions defined by fl3.1ip and fl3.12p . re¬ 
spectively. The spline Galerkin method fl3.ip is stable if and only if the operators 

R^ :RxR ^Rx R, 

R^ := 


/ I A(n,fi \ / 0 W M 0 \ (3.13) 

V AM I 0 )\ 0 M ) 

are invertible for all r G Air- 

Proof. By Corollary I3.1l the sequence (PnBrPn) is stable if and only if the coset 
(PnBrPn) + is invertible. Moreover, since Tp of fl2.ip is a compact operator, the 
sequences (P^ArPn) and {PnB-pPn) belong to the same coset (P^ArPn) + of the 
quotient algebra . However, by a version of the Allan’s Local Principle [T] 

for real C*-algebras [Til Theorem 1.9.5], the coset (P^ArPn) + is invertible if 
and only if for every r G P the coset {PfAY.^Pf) + is invertible in the corre¬ 
sponding algebra B^"" j . Therefore, the stability of our operator sequence will be 
established if we manage to show the invertibility of all cosets {PfAr.,.Pf) + 
r G P. Let us start with the case where r is not a corner point of P. If r ^ Air, 
then 6j. = tt, and straightforward calculations show that and Kj. are the zero 
operators. Hence, Ar is just the identity operator I in the corresponding space, so 
that PfArPf = Pf. The sequence {Pf) is obviously stable so that the corresponding 
coset (P^) + is invertible. 
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Consider next the case where r G Air- Note that by [71 Theorem 2.2] the 
operator A^ is invertible on the space L‘^{T). Therefore, by Corollary 13.21 the coset 
{P^ArPn) + is invertible in if and only if the sequence (P^ArP^) is 

stable. However, the stability of this sequence is equivalent to the stability of the 
sequence {'I>{P^At-P^)), where mapping is dehned by 03.51) . Consider also the 
operators A„ : —)■ P dehned by 

■^n ^ • • • ) )• 

By Lemma [2T] these operators are bounded and continuously invertible. Set A_„ := 
A“^ and note that the sequence (\1 /(P^At-P^)) is stable if and only if so is the 
sequence (P^), where 

Rl = diag(A„, A„) ■ d'(P,];A^P,];) ■ diag (A_„, A_„) :P xP x P. 


From the dehnition of the mappings d' and A±„ one obtains that the operators P^ 
have the form 


with the operators Ai^ ’ : P ^ P dehned according to the relations 

03.9p . However, these operators do not depend on the parameter n at all. Really, 
consider the matrix representations of the operators A^’^\ B^’^\ It 

follows from 03.10p that the entries aiq of the corresponding matrices are 

ds 


^pq — 


-f 

^ Jk+ 

-f 

^ Jr+ 


K(j)qn{a)(j)in{(j) da = 

u + I 


f—^ (pins — q)—(pina — 1) da 
\sJ s 


r-d+l 


= 


(p{t) 


t + q 

u + I 
t + q 

d+l 


, / du ,, . , 
(p{u) - (p{t) dt 


u + q' 


du 


'0 


{ndVn(p{u)) —^— {i^dVn(p{t)) dt 
u + q 

u ^, , du 

<p{u) -;- dt, 


t + q 


u + q 


hence these operators are independent of n. Moreover, = 0 and = 

^ 22 ’^^ = /. Combining all the above representations, one obtains that the operators 
P^ do not depend on the parameter n. Therefore, (P^) is a constant sequence and 
it is stable if and only if any its member, say R\, is invertible. It remains to observe 
that P"^ = P[, which completes the proof. ■ 
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4 Numerical approach to the invertibility of local 
operators 


As was already mentioned, there is no efficient analytic method to verify the 
invertibility of the local operators R'^. On the other hand, numerical approaches 
turn out to be surprisingly fruitful. Recall that the operators i?’", r G Adr do 
not depend on the shape of the contour F but only on the relevant angles 6^ and 
fir- Therefore, for contours having only one corner point. Theorem 13.21 can be 
reformulated as follows. 

Corollary 4.1 If t is the only corner point of the contour T, then the operator R'^ 
is invertible if and only if the Galerkin method {PnBrPn) is stable. 

Thus in order to determine the critical angles, i.e. the opening angles 6 for which 
the operators R'^ are not invertible, one can consider the behaviour of the spline 
Galerkin methods on special contours. A family of such contours Tf, 6 G (0, 27r), 

T^ := {f G C : f = 7i(s) = sin(7rs) exp(*6'(s — 0.5)), s G [0,1]} 

has been used in to study the local operators of the Nystrom method for 

Sherman-Lauricella and Muskhelishvili equations. Changing the parameter 6 in 
the interval (0, 27r), one obtains contours located at the origin and having only one 
corner of various magnitude. In the present paper, we use the same contours to 
detect the critical angles of the spline Galerkin methods. It is worth mentioning 
that the operator R'^ depends not only on 9r but also on the angle jSr between the 
right semi-tangent to the contour Tf at the point r and the real line M. However, 
numerical experiments conducted for both the Nystrom and spline Galerkin methods 
show that, in fact, the angle (dr does not influence the invertibility of the operator R'^. 
This opens a way for verifying the results obtained for contour T® by conducting 
similar tests for equations on contours with two or more corner, all of the same 
magnitude. To this end, we will use another contour T®, which is the union of two 
circular arcs with the parametrization 

7 i(s) = —0.5 cot(O.50) -|- 0.5/sin(0.56') exp( 26 *(s — 0.5)), 0 < s < 1, 

72 ( 5 )= 0.5cot(O.50) — 0.5/sin(0.56>) exp(f0(s — 0.5)), 0 < s < 1. 

To hnd the angles of instability, the interval [O.lvr, l.Ovr] has been divided by the 
points 6k '■= 7r(0.1-1-0.01/c) and for each opening angle 6k we constructed the matrices 
of the corresponding approximation operators for the Galerkin methods based on 
the splines of degree d = 0, d = 1 and d = 2. Note that we consider Galerkin 
methods for two choices of u, namely for n = 128 and n = 256, and the integrals 
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Figure 5: Condition numbers vs. opening angles in case n = 128. From row 1 to row 
3: splines of degree 0, 1 and 2, respectively. Left column: one-corner geometry, right 
column: two-corner geometry. 


arising in the equation fl2.2p and in the method 02.61) have been approximated by 
quadrature formulas 02.9p . O2.10p . Further, to verify the stability of the method, 
for each angle 6k we compute the condition numbers of the corresponding matrices 
and the results of these computations are presented in Figures [MU where possible 
presence of peaks might indicate critical angles. Thus it seems that inside of the 
interval (O.Itt, I.Qtt) neither of the Galerkin methods based on splines of degree 0,1 
or 2 has critical angles. This differs from the Nystrom method, where critical angles 
have been discovered for both Sherman-Lauricella and Muskhelishvili equations [SI 
[9]. Contrariwise, information about the critical angles at the interval ends is not so 
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conclusive. Thus in the case n = 256, the computation of the condition numbers for 
both one and two corner geometry shows that for the Galerkin method based on the 
splines of degree zero there can be a critical angle at the right end of the interval 
mentioned. 

For splines of the degree d = 0 and d = 1, the one and two corner geometries 
give contradictory results (see Figure E]). To clarify the situation one has to rehne 



Figure 6: Condition numbers vs. opening angles in case n = 256. From row 1 to row 
3: splines of degree 0, 1 and 2, respectively. Left column: one-corner geometry, right 
column: two-corner geometry. 

the mesh {6k} and essentially increase the dimension of the matrices used. Note 
that while discovering a suspicious critical angle for n = 256, we rehned the mesh 
{9k} in a neighbourhood of that angle by reducing its step to O.OOlvr, and calculated 
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the condition numbers for the corresponding Galerkin methods with n changed to 
512. This allows us to show that, in fact, there are no critical angles in the interval 
mentioned. However, the computing time increases drastically. 

The numerical experiments are performed in MATLAB environment (version 
7.9.0) and executed on an Acer Veriton M680 workstation equipped with a Intel 
Core i7 vPro 870 Processor and 8GB of RAM, and it took from one to two weeks of 
computer work in order to obtain every single graph presented in Figure |5l [6]or[3 


16r 

14 

12 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

0.1lt<e<1.9l[ 


16r 

14 

12 



qI-1-1-1-1-1-1-1-1-1-1 

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

0.lJ!<e<1.9Jl 






Figure 7: Condition numbers vs. opening angles in case n = 256 and n = 512 in 
neighbourhoods of suspicious points. From row 1 to row 3: splines of degree 0, 1 and 
2, respectively. Left column: one-corner geometry, right column: two-corner geometry. 
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